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ABSTRACT 

We show how future measurements of the SZ effect (SZE) can be used to constrain the 
cosmological parameters. We combine the SZ information expected from the Planck 
full-sky survey, N{S), where no redshift information is included, with the N{z) ob- 
tained from an optically-identified SZ-selected survey covering less than 1 % of the 
sky. We demonstrate how with a small subsample (w 300 clusters) of the whole SZ 
catalogue observed optically it is possible to drastically reduce the degeneracy among 
the cosmological parameters. We have studied the requirements for performing the 
optical follow-up and we show the feasibility of such a project. 

Finally we have compared the cluster expectations for Planck with those expected for 
Newton- XMM during their lifetimes. It is shown that, due to its larger sky coverage, 
Planck will detect a factor ^ 5 times more clusters than Newton- XMM and also with 
a larger redshift coverage. 
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1 INTRODUCTION 



Clusters of galaxies have been widely used as cosmological 
probes. Their modeling can be easily understood as they 
are the final stage of the linearly evolved primordial density 
fluctuations (Press & Schechter 1974, hereafter PS). As a 
consequence, it is possible to describe, as a function of the 
cosmological model, the distribution of clusters and their 
evolution, the mass function, which is usually used as a cos- 
mological test (Bahcall & Cen 1993, Carlberg et al. 1997, 
Bahcal & Fan 1998, Girardi et al. 1998, Rahman & Shan- 
darin 2000, Diego et al. 2000) Therefore, a detailed study of 
the cluster mass function will provide us with useful infor- 
mation about the underlying cosmology. 
Following this idea, several groups have tried to constrain 
the cosmology by using the information contained in the 
mass function. They compare the observational mass func- 
tion with the theoretical one given by the PS formalism 
(Carlberg et al. 1997, Bahcall & Fan 1998) or with simulated 
ones from N-body simulations (Bahcall & Cen 1993, Bode 
et al. 2000). The method has shown to be very useful but it 
is limitted by the quality of the data. Cluster masses are not 
very well determined for intermediate-high redshift clusters 
and even for low redshift ones the error bars are still sig- 
nificant. The standard methods to determine cluster masses 
(velocity dispersions, cluster richness, lensing, X-ray surface 
brightness deprojection) usually give different answers. It is 



believed that the best mass estimator is the one based on 
lensing (e.g. Wu & Fang 1997, Allen 1998) but the number 
of clusters with masses measured with this technique is too 
low to make this method a reliable technique to build a sta- 
tistically complete mass function although several attempts 
have been made (Dahle 2000). 

Instead of the mass function, it is possible to study the clus- 
ter population through other functions like the X-ray flux 
or luminosity functions or the temperature function. The 
advantage of these functions compared with the mass func- 
tion is that in these cases the estimation of the X-ray fluxes, 
luminosities, or temperatures of the clusters is less affected 
by systematics than the mass estimation from optical data. 
The drawback, however, is that to build these functions from 
the mass function, a relation between the mass and the X- 
ray luminosity, or flux or temperature is needed. These rela- 
tions are known to suffer from important uncertainties which 
should be taken into account. These uncertainties have their 
origin in the intrinsic scatter of these relations as well as in 
the quality of the observational data used to build them. 
There are three basic wavebands in which galaxy clusters 
are observed; optical (galaxies bounded to the cluster), X- 
ray (bremstrahlung emission from the very hot intra-cluster 
gas), and more recently in the mm waveband (Sunyaev- 
Zeldovich effect, SZE). 

The flrst clusters were observed with optical telescopes and 
also the flrst cluster catalogues were based on optical obser- 
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vations (Abell 1958; Abell, Corwin & Olowin (AGO) 1989; 
Lumsden et al. 1992 (EDCC); Postman et al. 1996 (PDCS); 
Carlberg et al. 1997b (CNOC)). However the optical identi- 
fication of a cluster is not a trivial task. First, it is not easy 
to define the cluster limits or cluster size. When observed 
in the optical band, a cluster appears as a group of galaxies 
which are bounded by the common gravitational potential. 
However, in the outer parts of the cluster there can be some 
galaxies for which it is not clear to what extent they are 
bounded to the cluster or, on the contrary, they are field 
galaxies. 

Optical identification of clusters has other more important 
shortcomings. They specially suffer from projection effects, 
that is, there can be some galaxies in the line of sight which 
are not bounded by any gravitational force but they can 
appear as a bounded system because their images are pro- 
jected onto a small circle centered on the cluster position. 
The best way to reduce this effect is by means of spectral 
identification. However, these kind of observations are time 
consuming and, when observing distant clusters, spectral 
identification is only feasible for the most luminous galaxies 
in the cluster. 

These problems are reduced when the cluster is observed 
in the X-ray band. In this band clusters appear as lumi- 
nous sources due basically to the continuos bremstrahlung 
emission coming from the very hot (T ~ 10*^ K) intracluster 
gas. The same intracluster emission can be used to deter- 
mine the cluster gas size. Moreover, projection effects are 
weaker in this case, thus X-ray surveys are very efficient in 
detecting clusters. However, with X-ray surveys it is diffi- 
cult to detect clusters further than z « 1. There are two 
reasons for this. One lies in the fact that the X-ray cluster 
emissivity is proportional to where Ue is the electron den- 
sity, that is, the emissivity drops very fast from the center 
of the cluster to the border and only the most dense cen- 
tral parts of the cluster will originate a substantial X-ray 
emission. This makes very difhcult the detection of distant 
clusters for which the X-ray emission is concentrated in the 
central parts of the cluster and therefore the apparent angu- 
lar size will be very small. Moreover the X-ray flux declines 
as D^^ = (_Dm(l -I- z))~^ {Dm is the comoving distance and 
Dl the luminosity distance). This selection function limits 
the redshift at which one cluster can be observed by actual 
X-ray detectors which are blind to the earlier stages of clus- 
ter formation (2: <: 2). 

In the mm waveband the situation is quite different. The 
SZE surface brightness goes as n^. Therefore, when observ- 
ing the clusters at mm wavelengths, it would be easier to 
detect the most distant ones because they will show a larger 
apparent angular size since the brightness profile drops more 
slowly than in the X-ray case. In X-rays, the total flux de- 
cays as D]^^ ■ Meanwhile, the integrated SZE flux goes as 
Da = {Dm/{1 + z))~ (Da being the angular diameter dis- 
tance). Da grows more slowly than Dl and even decreases 
after a certain redshift 2: « 1. Therefore, the SZE flux drops 
more slowly with distance (or even increases) in the case of 
the SZE than in the X-rays case. Another advantage is that, 
as in the X-ray situation, galaxy clusters observed through 
the SZE does not suffer (almost) from projection effects. All 
these reasons make the SZE the preferable way of observing 



distant clusters. 

Our interest on the SZE is twofold. First, it can be consid- 
ered as a contaminant of the cosmological signal (CMB) and 
therefore a good knowledge of this effect is required in order 
to perform an appropriate analysis of the CMB data. But 
it can also be considered as a very sensitive tool to measure 
the mass-space cluster distribution. In this paper we will 
concentrate our effort in this second aspect. 
Planned CMB surveys will also be sensitive to the SZE dis- 
tortion induced by galaxy clusters (MAP, Planck). These 
surveys will cover a wide area of the sky and they are ex- 
pected to detect the SZE signature for thousands of clusters. 
Furthermore, proposed and undergoing mm experiments will 
measure the SZE for hundreds of clusters in the near future 
(e.g. AmiBa, Lo et al. 2000; LMT-BOLOCAM, Lopez-Cruz 
& Gaztaiiaga 2000; CBI, Udomprasert et al. 2000; AMI, 
Kneissl et al. 2001; or the interferometer proposed in Holder 
et al. 2000) The cosmological possibilities of these new data 
sets are very relevant as they will have a statistically large 
number of clusters and they will improve the redshift cov- 
erage significantly. For a realistic prediction of the power of 
future SZE surveys we should consider the detector charac- 
teristics of one of these planned experiments. In this work 
we will focus our attention on the expected SZE detections 
for the Planck mission and we will study the possibilities 
of these data to probe the cosmological model. Planck will 
observe the whole sky at 9 mm frequencies (including those 
where the SZE is more relevant) and with angular resolu- 
tions ranging from 5 arcmin to 33 arcmin FWHM (see fig. ^ 
below). The observation of the SZE at different frequencies 
will make easier the task of identifying clusters because of 
the peculiar spectral behaviour of the SZE signal which can 
be very well recognized with the nine Planck frequencies. As 
can be seen from figure (|lj), the best channels are the ones 
at 100 GHz {x = 1.76), 143 GHz {x = 2.5) and 353 GHz 
{x = 6.2) together with the channel at 217 GHz {x = 3.8) 
where the thermal SZE vanishes. A cross correlation of these 
channels, including the knowledge of the spectral shape fac- 
tor, will allow to discriminate among the SZE, foregrounds 
and CMB. 

As we will show in the next sections, the cluster population 
at high redshift is much more sensitive to the cosmological 
parameters than the cluster population in the local Universe 
{z ^ 0.5). Thus, the information at high redshift is crucial 
to determine the underlying cosmology and the SZE can be 
the tool to get that information. 

The structure of this paper is the following; in section ^ 
we recall some of the basics of the SZE and some useful 
definitions which will be of interest in subsequent sections. In 
section ^, we show how through the SZE we can investigate 
the cluster population and evolution. We also show in section 
how future SZE surveys should be complemented with 
optical observations of a small subsample of SZE-selected 
clusters in order to provide redshift information of those 
clusters; in this way, it is possible to reduce the degeneracy in 
the cosmological models describing the data. We apply this 
idea in section]^ to simulated Planck SZ data. In that section 
we obtain some interesting results about the possibilities 
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of those future data. Finally we discuss our results and we 
conclude in section ^. 



2 THE SUNYAEV-ZEL'DOVICH EFFECT 

Since Sunyaev & Zel'dovich predicted that clusters would 
distort the spectrum of CMB photons when they cross a 
galaxy cluster (Sunyaev & Zel'dovich 1972), several detec- 
tions of this effect have been made. At present, the number 
of clusters with measured SZE is small because they are 
limitted by the detector sensitivity (a typical SZE signal is 
of the order of 10"" in AT/T with T ^ 2.73 K). However, 
the sensitivity in the detectors is getting better and better 
and in a few years the number of detected clusters through 
the SZE will make this effect an important tool to explore 
clusters at high redshift. 

When CMB photons cross a cluster of galaxies, the spectrum 
suffers a distortion due to inverse Compton scattering. The 
net distortion in a given direction can be quantified by the 
cluster Comptonization parameter, yc, which is defined as 



yc 



Tne{l)dL 



(1) 



The integral is performed along the line of sight through the 
cluster. T and are the intracluster electron temperature 
and density respectively; or is the Thomson cross section 
which is the appropriate cross section in this energy regime 
and fcs is the Boltzmann constant. 

The flux and the temperature distortion are more widely 
used than the Compton parameter because these are the 
quantities which are directly determined in any experiment 
when observing clusters at mm frequencies. The distortion 
in the background temperature is given by: 



(2) 



5T_ 
~T 

where x = hv/kTcMB is the frequency in dimensionless 
units and g{x) = xcoth{x/2) — A is the spectral shape factor. 
There is a second kind of distortion due to the bulk velocity 
of the cluster. It is known as the kinetic SZE. However, this 
contribution is much smaller than the thermal contribution 
and we will not consider that in this work. 
As can be seen from the previous equation, there is no red- 
shift dependence on the temperature distortion, i.e, the same 
cluster will induce the same distortion in the CMB temper- 
ature, independently of the cluster distance (except for rel- 
ativistic corrections). The only redshift dependence of the 
total SZE flux is due to the fact that the apparent size of 
the cluster changes with redshift. 
The total SZ flux is given by the integral ; 



Stotal{v) 



AI{iy,e)dQ, 



(3) 



where the integral is performed over the solid angle sub- 
tended by the cluster. AI{v,d) is the change in intensity 
induced by the SZE and is given by AI{i',9) — Iof{x)yc 
where Jo = {2h/c^){kBTcMB /hf and 



a:(e" + 1) 
(e^ - 1) 



(4) 



AI = I„f(x)y„ 1. . (zh/c')(k,T„,/h)- 



Planck-SZE 




Figure 1. Spectral shape of the thermal SZE. Planck channels 
are between 30 GHz and 800 GHz (black dots), including also 
the 217 GHz channel where the thermal SZE vanishes. The best 
channels to detect the SZE decrement/increment will be those of 
the HFI instrument at 100 GHz {x = 1.7), 143 GHz {x = 2.5) 
and 353 GHz {x = 6.2) respectively together with the one at 217 
GHz (x = 3.8). At 100 GHz (x = 1.7) there are two channels: the 
LFI (solid line) and the HFI (dotted line). 



is the spectral shape (see fig.^). 

In fig. hi we show the spectral shape of the SZE f{x). In the 
same plot the 9 frequency channels from Planck are shown 
as black dots. The amplitude of the lines are proportional to 
the SNR per resolution element in each channel (assuming 
that the clusters are unresolved). The SNR depends on the 
spectral shape factor f{x) since the signal is proportional 
to f{x) (see eq. ^ below) and the sensitivity per resolution 
element (cr); SNR — S/N oc f{x)/a which is proportional 
to the amplitudes represented in fig. ^. 

Eq. can be easily integrated assuming that the clus- 
ter is isothermal. Then the integral can be reduced to 
oc J dfl J ne{9)d£ which can also be transformed into 

M„ 



D-\z)J^M. M, h 



D-\z) J dVn^e) = D-2(z) 
and nip are the total mass, baryon fraction and the proton 
mass respectively. In this approach, we have assumed that 
the gas is only composed of ionized Hydrogen. Finally we 
get: 



SsZE = 



3.781/(a)/ftrMi5 
DUz) 



(5) 



In the previous expression T is given in Kelvin, M15 in 
10^^ Mq and Da (z) in Mpc. In these units, the flux is given in 
mjy (1 mjy = 1.0 x 10"^'" erg s~^ Hz~^ cm~^). If we now 
introduce the h dependency in M15 (10^^ Mq), Da{z) 
{h~^ Mpc) and fb ~ 0.06/i~^ then the flux is given also in 
mJy units. 

Eq. (j^ tell us that the total flux is given basically by three 
terms: temperature, mass and angular-diameter distance to 
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the cluster. As a difference with the X-ray flux, where in 
the integral over the cluster volume one should consider n1 
instead of just rie , the total SZE flux does not depend on the 
cluster profile. Therefore, for an isothermal cluster, no as- 
sumption about the electron density profile is needed when 
computing the total SZE flux. This aspect is very important 
since for the Planck resolution, only a small percentage of 
the clusters (« 1 %) will be resolved and as a first approxi- 
mation we can consider that all the fiuxes can be computed 
from eq. (^. This fact will simplify significantly our calcula- 
tions and simultaneously will reduce the uncertainty due to 
the lack of a precise knowledge of the cluster density profile. 
The second important consequence that we can remark from 
the previous equation is about the angular-diameter distance 
dependence of the total flux. In X-rays, the total flux de- 
pends on D2^{z). If we compute the ratio D1^{z)/D\{z) it 
goes like (1-1- z)"^, that is, the Ssze drops a factor (1-f 
slower than the Sx ■ Thus for high redshift clusters it is evi- 
dent the advantage of using mm surveys (SZE) as compared 
with the X-ray ones. 

If the cluster is unresolved by the antenna, then the cluster 
total flux (eq. |^ can be considered contained in the resolu- 
tion element of our antenna. If on the contrary, the cluster 
is resolved, then the previous approximation is not adequate 
anymore and we should integrate the surface brightness over 
the cluster solid angle. 

The gas density proflle can be well fitted by a /3-model. 

n,{r)^no{l + {r/r,f)-^^'\ (6) 

where Uo is the central electron density (rio ~ 2.0 x 10~^(1 + 
zYh^^'^e~ /cm^) and rc is the core radius of the cluster. The 
observed values of /3, obtained from X-ray surface brightness 
profiles, typically range from 0.5 to 0.7 (Jones & Forman 
1984, Markevitch et al. 1998). 



3 SZE AS A PROBE OF THE 
COSMOLOGICAL MODEL 

As we mentioned in the introduction, an estimate of the clus- 
ter mass function can be used to constrain the cosmological 
parameters. However such estimate is affected by the sys- 
tematics in the mass estimators. It is possible to use other 
functions to explore the cluster population like for instance 
the temperature function cP'N/dTdz (Henry & Arnaud 1991, 
Eke et al. 1998, Viana & Liddle 1999, Donahue & Voit 1999, 
Blanchard et al. 2000, Henry 2000). The connection between 
the mass function and the temperature function is the T—M 
relation, 

(fN _ (fN dM 

dTdz ^ dMdz~dT ^ ' 

Temperature estimates are more reliable than mass esti- 
mates. Therefore the temperature function should be less 
affected by the systematics than the mass function. How- 
ever, the scatter in the T — M relation introduces new un- 
certainties which should be taken into account. Similar prob- 
lems have the X-ray luminosity function and the X-ray fiux 
function which have been widely used in the literature to 
constrain the cosmological parameters (Mathiesen & Evrard 
1998, Borgani et al. 1999, Diego et al. 2000). In those cases 



Table 1. Best ACDM (fiat) and OCDM (A = 0) models. Both 
models are undistinguishable in the redshift regime of the data 
considered in Diego et al. (2000). 
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new uncertainties appear due to the scatter in the Lx — T 
relation. In a previous work (Diego et al. 2000), we have con- 
sidered these ideas and applied them to constrain the cosmo- 
logical parameters by fitting the model to different data sets: 
mass function, temperature function. X-ray luminosity and 
flux functions. That work was innovative in the sense that 
we considered all the previous data sets simultaneously in 
our fit (and not just one as usual) and we looked for the best 
fitting model to all the data sets considered. Also important 
is to remark that in that work we considered the cluster scal- 
ing relations as free-parameter relations. These two points 
are very relevant because when fltting cluster data sets it is 
important to check that the best fitting models agree well 
with other data sets. Otherwise, we should reexamine the 
assumptions made in the model. One of these assumptions 
is, for instance, the cluster scaling relations T — M or Lx — T. 
These relations are commonly fixed a priori. However they 
are known to suffer from important scatter as we have shown 
in that work. In fact, we found that not all the scaling rela- 
tions considered previously in other works were appropriate 
to describe conveniently several cluster data sets in a simul- 
taneous fit. By fitting our free-parameter model to all the 
data sets, we obtained, not only the cosmological parame- 
ters, but also the best scaling parameters. In that work we 
found as our best fitting parameters the ones given in ta- 
ble 0. The first 3 columns correspond to the cosmological 
parameters and the others are for the T — M relation; 

T,as=ToMWi- + z)^, (8) 

where M15 is the cluster mass in 10^^/i~^Mq and To, a, 
and tp were free parameters. We found that the two mod- 
els in table ^ were undistinguishable when they were com- 
pared with the data. In fact, both models agree well with the 
data considered and there was no way to determine which 
one was the best. In order to distinguish them, more and 
better quality data are needed. Undergoing X-ray experi- 
ments (CHANDRA, Newton-XMM) will help to do that as 
well as current and proposed optical surveys (e.g. Gladders 
et al. 2000). The models which were undistinguishable in 
the low redshift interval will show different behaviour at 
higher redshift. However, as we mentioned in the introduc- 
tion, the waveband in which distant clusters will be detected 
the most, is not the optical nor the X-ray band, but the mm 
band. The usefulness of the SZE as a tool to measure cosmo- 
logical parameters motivated many works in the past (e.g. 
Silk & White 1978, Markevitch et al. 1994, Bartlett & Silk 
1994, Barbosa et al. 1996, Aghanim et al. 1997) as well as 
in the present time (Holder et al. 2000, Grego et al. 2000, 
Mason et al. 2001, this work). 



© 0000 RAS, MNRAS 000, 



SZE & Cosmology 5 



SZE (353 GHz) XCDM OCDM 

10' F ^ I II ^ ^ ^ ^ ^ 11 ^ ^ 1 



10° r 




Figure 2. The SZE integrated number counts of the cluster 
population predicted at mm wavelengths (353 GHz) for the best 
flat (solid) and open (dotted) models in table 1. In the plot three 
redshift shells are represented: top z < 1, middle z S [1, 2] and 
bottom 2 > 2. 

Using SZ data it is possible to estimate the cosmological pa- 
rameters by looking at some SZ derived function related with 
the mass function. We just need a measurement of a cluster 
quantity related with its mass. As we have shown in eq. 
the total SZ flux can be such a quantity. From that equa- 
tion it is possible to build a mm flux function (dN/dSmm) 
from a mass function (dN/dM). An interesting alternative 
to this approach can be found in Xue & Wu (2001). In that 
paper the authors have connected the mm flux function with 
the observed X-ray luminosity function through the Lx — T 
and T — M relations. If we consider a future mm experiment 
(like Planck) where thousands of clusters are expected to be 
detected through the SZE, probably a fit to the mm fiux 
function could be able to distinguish the models in table |l| 
In fig. 1^ we show the expected number of clusters with total 
flux above a given flux for the two previous models. These 
are integrated fluxes, i.e we assume that the clusters are not 
resolved and that their total flux falls into the antenna beam 
size. This assumption is appropriate when the antenna size 
is above several arcmin as it is the case for the Planck satel- 
lite where its best resolution is FWHM = 5 arcmin. 
In fig. ^ we did not show the whole (0 < z < oo) inte- 
grated N{> S) curve since it looks quite similar to the two 
curves on the top of the figure and they still look undistin- 
guishable. To compare them in a more realistic way we have 
performed a Kolgomorov-Smirnov test where we have com- 
pared the A'^(> S) curve for a realization (over 2/3 of the 
sky) of the ACDM model with the mean expected values of 
both OCDM and ACDM models. The KS test was unable 
to distinguish between the models. As can be seen from the 
figure, both models predict the same number of clusters in 
the z < \ redshift interval (top). That situation was similar 
in Diego et al. (2000) where most of the cluster data were at 
low z and it was not possible to discriminate between both 
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Z 

Figure 3. The SZE differential number counts of the cluster 
population predicted for the Planck mission at mm wavelengths 
(353 GHz) for the best flat (solid) and open (dotted) models in 
table 1^ as a function of redshift. 

models. 

But the situation changes at redshifts z > 1. In the red- 
shift bin z G (1,2) (middle) the differences between both 
models are significant. These differences are increased when 
we compare the models in the redshift bin z > 2 (bottom) 
where in the OCDM model two orders of magnitude more 
clusters are expected than in the ACDM case above 5* = 30 
mjy. This flux [S = 30 mjy at 353 GHz) is expected to be 
the flux limit of Planck (see Appendix A where we estimate 
that limit based on MEM residuals). However, this limiting 
flux will depend on the method used to identify such clus- 
ters in the maps. The final method to be used with Planck 
is still under development. MEM methods work very well 
(Hobson et al. 1998, Hobson et al. 1999) being at present 
the preferred one. 

In fig. P we show the number of clusters with fluxes above 
that limiting flux of Planck as a function of redshift. Again 
the differences are not significant below z ^ 1 but they are 
quite relevant above that redshift. By looking at figs. (^ 
(m), we see that the differences in the cosmological mod- 
els are more evident at high redshift. In order to discrim- 
inate among the cosmological models one should consider 
not only the cluster population at low redshift (normaliza- 
tion) but the cluster population at high redshift (evolution) 
as well. A recent application of this idea can be found in 
Fan & Chiueh (2000) where the authors study the function 
r = N{z < 0.5) /N{z > 1) as a function of the limiting flux. 
That work suggests that a limitted knowledge of the redshift 
of the cluster would allow to constrain the cosmological pa- 
rameters. 

Both fig. ^ and fig. |^ suggest that with future SZE data it 
will be possible to go further on the determination of the 
cosmological parameters. From all these plots, it is evident 
that Planck (together with redshift information for a small 



© 0000 RAS, MNRAS 000, 



6 J.M Diego et al. 



subsample, see next section) would be able to discriminate 
between the models which previously were undistinguishable 
when they were compared with present X-ray and optical 
data. Furthermore, the accuracy in the free parameters ob- 
tained in previous works will be increased with these new 
data. 



3.1 Montecarlo simulations: Planck vs 
Newton-XMM 

The previous discussion was based on theoretical expecta- 
tions of the mean number of SZE detections expected at dif- 
ferent redshifts. We want to go further by computing Mon- 
tecarlo (Press-Schechter) realizations of the expected SZE 
on a specific patch of the sky. In order to compare the mm 
and X-ray bands we will also compare the expected SZ de- 
tections in this area for Planck with those based on X-ray 
expectations for Newton-XMM in the same sky patch. 

The simulations were done over a patch of the sky of 
12.8° X 12.8° and with a pixel size of 1.5' x 1.5' fihered with 
a FWHM of 5 arcmin and at the frequency of 353 GHz, fol- 
lowing the characteristics of the 353 GHz channel of Planck. 
The parameters of this simulation corresponds to the ACDM 
model of table 1. 

The total number of clusters is about ~ 20000 in one of these 
simulations. The mean Comptonization parameter is well 
below the FIRAS limit (~ 10"^ compared with 2.5 x 10"^, 
Fixsen et al. 1996). The resulting distribution of clusters is 
shown in figure ^. The allowed range of simulated masses was 
3.0 x W^h-^MQ < M < 1.0 X The lower limit 

is at the frontier between a small cluster and a galaxy group 
and it is the mass corresponding to a cluster with temper- 
ature T ~ 1 keV which can be considered as the minimum 
temperature for a virialized cluster. 

In fig. ^ each point corresponds to one cluster with mass Al 
and redshift z. This distribution is in agreement with the 
observational constraint given by the detection of at least 
three clusters with masses above 0.8 x IO^^Mq and z > 0.5 
(Bahcall et al. 1998). 

The most massive cluster is at redshift 0.66 with a mass 
M = 1.09 X W^^h-^MQ ahhough this was an unusual real- 
ization. In a normal one of this size, the most massive cluster 
is usually well below z ~ 0.5. Clusters marked with an open 
circle correspond to those with a total fiux above 30 mjy and 
according to our criterion, these clusters would be detected 
by Planck (see Appendix A). There are 185 clusters above 
this limit and only one of them would be observed above 
z = 1 in this simulation. As a comparison we show the same 
picture but corresponding to the OCDM model in table 1. 
(fig. 1^. In this case ~ 15 clusters are detected above redshift 
z = 1. This comparison demonstrates how a small region of 
the sky can show up the differences between the two models 
in table 1. Both models are nearly undistinguishable below 
z ~ 0.7 but they difi^er significantly above that redshift. We 
will come later to this point in section |^. 

As a comparison with the number of detected clusters ex- 
pected with Planck, we show in figs. §and| the clusters 
expected to be detected by Newton-XMM (big solid circles, 
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Figure 4. ACDM model. Each cluster is represented as a dot in 
this plot. Clusters marked with a big black dot have fluxes (in the 
0.5 — 2 keV band) above the limiting flux Sx = 1-4 X 10~^* erg 
cm~^ and they are the clusters expected to be detected by 
Newton-XMM in the XCS (Romer et al. 1999). Cluster marked 
with an open circle are those with a SZE flux above 30 mJy at 
353 GHz and they are the clusters expected to be detected by 
Planck. Solid lines represent the masses for which the flux equals 
the limiting flux of Planck (30 mJy), and the XCS. Dotted lines 
show the maximum cluster mass expected in a solid angle equal 
to 2/3 of the sky (top) and (12.8°)^) (bottom). Finally the heavy 
solid line indicates the selection function for an hypothetical op- 
tical survey made with a 10-m class telescope and t wo hours of 
observation (8 bands and 900s per band. See section 



Sx > 1.5 X 10"" erg s"^ cm'^ in the 0.5 - 2 keV band, 
Romer et al. 1999). 

Solid lines represent the minimum mass as a function of red- 
shift for which the corresponding fiux is above the limiting 
fiux in Planck and Newton-XMM respectively, that is, they 
are the selection functions for both missions. 
Although Newton-XMM will see many more clusters at low 
redshift, however Planck will be more sensitive to those clus- 
ters at high redshift. The reason is that the corresponding 
selection functions are different in both cases. The X-ray fiux 
goes as Sx oc D^'^{z) and to detect clusters deeper in red- 
shift they must be more luminous (more massive) in order 
to compensate the monotonous increasing function Dl{z). 
On the contrary, the SZE fiux goes like Ssze oc Da'^{z) 
and at redshift ~ 1, the angular diameter distance reaches a 
maximum and after that it starts to drop to smaller values. 
Therefore the masses needed to provide a particular fiux, 
Ssze = 30 mJy, can be smaller at redshift z > 1 than the 
required masses at z « 1. 

The dotted lines in figs. (^) and (^) indicate the maximum 
expected mass as a function of redshift in the simulation for 
the specific solid angle 2/3 of the sky (top) and 12.8^ deg^ 
(bottom). The point where the selection function cross these 
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Figure 5. OCDM case. Note the differences in the cluster dis- 
tribution above the coasting phase redshift Zc ~ 0.7 (see text). 

maximum mass lines corresponds to the maximum redshift 
expected in the sample. In both cases, Planck goes deeper 
than Newton-XMM. The difference can be even more em- 
phasized if the limiting flux for Planck is below the one con- 
sidered in this paper which was a very conservative one (see 
Appendix A). 

As noted in Romer et al. (1999), the XMM Cluster Survey 
(XCS) will cover ~ 800 square degree and will contain more 
than 8000 clusters. On the contrary Planck will observe the 
full sky, that is the sky coverage will be ~ 50 times larger 
than that of the XCS. In order to compare Newton-XMM 
with Planck this difference in the sky coverage should be 
considered. In fig. ^ we show the number of detected clus- 
ters expected by Planck and Newton-XMM (XCS) where 
we have taken into account the sky coverages and limiting 
fluxes expected for both missions. 

Now, the differences between Planck and Newton-XMM are 
evident. The large sky coverage of Planck together with the 
constancy of the SZE surface brightness with redshift, will 
allow this satellite to detect a much more significant number 
of clusters than Newton-XMM in the XCS. Also, we can 
conclude from this plot that, for the best fitting models in 
table 1, no clusters are expected to be detected above z ~ 3, 
neither for the OCDM nor the ACDM. However the cluster 
abundance between 2 = 1 and z = 2 will provide a definitive 
probe of the cosmological parameters and the cluster scaling 
relations. 

An important consequence of the previous plots (figs, ^and 
is that only a small portion of the whole sky would be 
needed to distinguish between the two models ACDM and 
the OCDM. This is an important point because only spectral 
identification of ~ 200 (not resolved) random selected clus- 
ters from the whole catalogue would be needed. However, it 
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Figure 6. Expectations for Planck (solid line OCDM, dotted 
lino ACDM) and the XCS (dot-dashed line OCDM, dashed line 
ACDM). 

is important to answer the question, what is the minimum 
number of clusters needed to distinguish the models in table 
1 at, for instance, the 3 cr level? We will try to answer this 
question in the next section. 



4 AN OPTICALLY-IDENTIFIED 
SZ-SELECTED CATALOGUE 

As we have seen in fig. ^ the information provided by an 
hypothetical N{S) curve (even if this curve corresponds to 
a full sky survey) will be insufficient to distinguish the two 
models considered in that figure. Redshift information will 
be needed in order to make the distinction. Different cos- 
mological models predict different cluster populations as a 
function of redshift. If we analyze the evolution of the cluster 
population with z, then we should be able to discriminate 
among those models. However, to study the evolution of the 
cluster population we need spectral identification of the clus- 
ters (or at least optical observations in several bands in order 
to get photometric redshifts) since the SZE does not provide 
any estimate of z. Performing these observations for an hy- 
pothetical full sky SZ catalogue would be a huge task but if 
only a small number of unresolved clusters need to be identi- 
fied then the work is significantly simplified. Now, we should 
ask the question, how small can our optically-identified sam- 
ple be if we want to distinguish between, for instance, the 
two models in table ^ To answer that question we have com- 
pared the curves N{Smm > 30 mjy, > z) for the two models 
in which we are interested. We require that at a given red- 
shift, both curves must be distinguished at a 3cr level, that 
is, we require the condition iV° - SVW = + 
(assuming Poissonian statistics), where A'^'^ and A'^^ are the 
number of detected clusters above a given z in the two cases 
OCDM and ACDM respectively. 

Since we know A'''"' and A'^ at each z, then we can com- 
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Figure 7. Required number of clusters to be observed in order 
to distinguish tlie OCDM and ACDM models. 

pute the required total number of clusters which should be 
observed in order to satisfy the previous condition at each 
redshift. fn fig |^ we show this calculation for three different 
selection criteria of the clusters. In each one of the lines we 
show the total number of clusters randomly selected from 
the catalogue (with the only condition that the total flux 
must be Smm > 100 mjy top, Smm > 30 mjy middle and 
30 mJy < Smm < 40 mJy bottom) which should be opti- 
cally observed in order to distinguish (at a 3(t level) A'''^ and 
A'^'^ at redshift z (i.e. the total number of observed clusters 
needed to have a 3a difference in A'^(> z) for the two mod- 
els). As can be seen from the plot, randomly selecting about 
300 clusters with Smm > 30 mJy from the full catalogue 
and determining the redshift for each one, will allow to dis- 
tinguish the two models at a 3a level above z « 0.6 just by 
looking at the A'^(> z = 0.6) curve. The explanation for this 
fact is given by the different evolution of the cluster forma- 
tion in both models. In the ACDM case there is a coasting 
phase (or inflection point in the acceleration parameter) at 
Zc ~ 0.7 that helps to form structure at this redshift. This 
phase is not present in the OCDM case where there is a red- 
shift Zc ~ 2 below which the collapse of linear fluctuations 
is inhibited by the fast expansion of the Universe. 
Choosing the selection criteria 30 mJy < Smm < 40 mJy, 
it is possible to reduce slightly the number of clusters to be 
identified. If on the contrary, only the brightest clusters with 
Smm. > 100 mJy are identified optically, then we would need 
a significantly larger number of clusters in order to make the 
distinction between the models. 

4.1 Cluster optical detection simulations 

Probably the most cost efficient way of identifying in the 
optical the galaxy clusters detected by Planck is using pho- 



tometric redshifts. Even a rough estimate of the photo-z 
allows to considerably reduce the background galaxy con- 
tamination and enhance the contrast density of the cluster. 
In addition, although the error in the photo-z of an individ- 
ual galaxy is usually ~ 0.1 at 2 < 1 (Bem'tez 2000), the total 
error in the cluster redshift will be 0.1 /y^Nci {Nd being the 
number of galaxies in the cluster). To estimate the feasibil- 
ity of detecting Planck SZ cluster candidates using optical 
imaging, we perform simulations based on empirical infor- 
mation. Since extensive data are only available for relatively 
low redshift clusters, this has the disadvantage of ignoring 
evolution effects. However, it has been shown that the evolu- 
tion of the cluster early type galaxies is not dramatic up to 
z m 1.3 (Rosati et al. 1999). Therefore, at worst, this makes 
the results obtained here a conservative, lower limit on the 
detectability of high redshift clusters, since any reasonable 
luminosity evolution would tend to make the cluster galax- 
ies bluer and brighter, increasing the number of detected 
galaxies with respect to the non-evolution case. 
Wilson et al. (1997) represent the luminosity function of 
bright galaxy clusters as the superposition of a Schechter 
function and a Gaussian for the brightest cluster galaxies. 
The spectral fraction for the cluster galaxies can be derived 
comparing the V — I colors of galaxies in A370 and CL 
1447-1-23 (Smafl et al. 1997) with those expected for El, Sbc, 
Scd and Im spectral templates (Coleman, Wu & Weedman 
1980). With the above luminosity function and spectral frac- 
tions (extended to 8 magnitudes fainter than M*) we gen- 
erate a bright galaxy cluster at z = 0.18. By redshifting 
this cluster, we generate a series of mock cluster catalogs 
containing the / band magnitude and the spectral type T 
from 2 = 0.2 to 2 = 2.0. The I— band magnitudes are 
transformed using k-corrections derived from the Coleman, 
Wu and Weedman templates. To model the surface number 
counts distribution of the clusters, we use a n{R) oc R~^'^ 
law, found by Vilchez et al. 2001 to agree well with the pro- 
jected galaxy density in the central regions of galaxy clus- 
ters. To link the optical results with X-ray and SZ quan- 
tities, we integrate the luminosity function of the cluster 
in the V— band, and assume M/L = 300, which leads to 
M ~ lO^'^M© for a A1689-fike cluster. 

We simulate the background galaxy distribution n{z, T, I) 
using the Hubble Deep Fields (Williams et al. 1996, Williams 
et al. 2000). For the redshifts, we use the spectroscopic re- 
suhs of (Cohen et al. 2000) for the HDFN and photo-2 
obtained using the BPZ code (Bemtez 2000) for the rest of 
the HDFN galaxies and those of the HDFS. 
Once we have I,z,T catalogs for all the galaxies contained 
in the field, we generate UBVRIJHK magnitudes using 
the above mentioned template library, enlarged for the HDF 
galaxies using two starsburst templates from Kinney et al. 
1996. Gaussian photometric errors are added to these ideal 
magnitudes using empirical relationships derived from real 
observations with 10m class telescopes, scaled by the square 
root of the exposure time needed to simulate a 900s per band 
observation. 

For cluster detection, we look for an overdensity of ellip- 
ticals with respect to the expected background population. 
The reason to do this, instead of using the whole cluster pop- 
ulation, is that the density contrast is much higher for this 
type of objects. Even a moderate cluster at 2 ~ 1 (Bem'tez 
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et al. 1999) conspicuously stands out against the relatively 
sparce numbers of field ellipticals (see also Gladders & Yee 
2000). Therefore, we estimate photometric redshifts for all 
the galaxies within an angular aperture corresponding to 
~ l/i~^Mpc, classify them into different spectral types, and 
construct a redshift histogram for the early types. The pres- 
ence of a 'spike' in the redshift histogram will indicate the 
existence of a cluster. The signal-to-noise of such detection 
is 

S/N^[N{z)-<N{z) >]/a,iz) 

where iJg{z) is the expected fluctuation in the galaxy num- 
bers within a redshift slice centered on z. Its value can be 
estimated as: 

ag{zf =< N{z) > +1^^M^ J w{ei2)dSidS2 

N{z) is the detected number of galaxies with an area S, 
< N{z) > is the expected average density and w(Si2) is 
the two-point correlation function within the redshift slice. 
For most of the redshifts considered, lMpc/i~^ corresponds 
to ~ 4 arcmin. Taking the amplitude of w{6i2) to be 
A « 7.6 X 10~^ within a. z = 0.2 slice (Brunner, Szalay 
& Connolly 2000), which is approximately the same redshift 
interval used here to detect the clusters, the value of cr{z) 
is roughly ag{zf ^< N{z) > (l + 3.27 x 10"^ < N{z) >). 
This number may be an understimate since the clustering 
strength of the early types is known to be higher than for the 
general galaxy population. The numbers below are based on 
a 3a detection limit as defined by the above equation. There 
are plenty of other methods, parametric and non-parametric, 
which will probably be more efficient in finding clusters (see 
e.g. Nichol et al. 2000). But again, we think that using a rel- 
ative simple procedure provides a good idea about the prac- 
ticality of this approach. A reasonable observing strategy 
will be to start with those clusters not detected in shallower 
imaging surveys, e.g. the SDSS catalog, and depending on 
the redshift/luminosity range of interest (e.g. low mass/low 
redshift clusters or high mass/high redshift ones), use only 
a small subsample of the BVRIJHK filter set mentioned 
above, which brackets the 4000A break at the required red- 
shift, and which would be enough to detect the early types. 
If one desires to reach a higher precision in the photo-z esti- 
mation, or wants to reach the limits shown in Figs ^ and ^ at 
all the redshift intervals, then the whole filter set should be 
used. The optical selection function in Figs ^ and |^ presents 
quite a jagged look. Apart from relatively smooth efi'ects 
like the cosmological dimming or the K-corrections which 
determine the general trend, the detectability of the clus- 
ters is significantly affected, at least at z < 1, by the rela- 
tive placement of the redshifted 4000A break with respect 
to the filter set, specially to the R and I band filters, the 
ones which go deeper for a fixed exposure time. When the 
break falls almost exactly in between these two filters, the 
photo-z precision is improved, whereas if the break is close 
to the central position of a filter the redshift error increases 
and the estimation is more easily affected by color/redshift 
degeneracies. Therefore, the relative exposure times and the 
filter choice should be optimized depending on the redshift 
interval that it is being targeted. 



5 ESTIMATING THE COSMOLOGICAL 
MODEL FROM N{S) AND N{Z) 

Previous works (Kitayama & Suto 1997, Eke et al. 1998, 
Mathiesen & Evrard 1998, Borgani et al. 1999, Diego et al. 
2000) have shown the power of X-ray surveys to constrain 
the cosmological model. From the previous sections we con- 
clude that also the SZE data can be used with the same pur- 
pose (see also Markevitch et al. 1994, Barbosa et al. 1996, 
Aghanim et al. 1997, Holder & Carlstrom 1999, Majumdar 

6 Subrahmanyan 2000, Fan & Chiueh 2000). 

As we have seen in the previous sections, both the N{S) and 
the N{> 30 mjy, 2) curves can be used to study the clus- 
ter population, being N{S) the curve having larger number 
of clusters (although with no z information) and A'^(> 30 
mjy, z) the curve which is more sensitive to the evolution of 
the cluster population. Following Diego et al. (2000) we will 
combine both curves to reduce the degeneracy. Some models 
which are compatible with the first curve will be incompat- 
ible with the second one and vice versa. Thus, those models 
will be rejected. 

Since, this kind of data is not available yet, we will check 
the method with two simulated data sets following the char- 
acteristics of the Planck satellite (section ^) for N{S). For 
the second curve we will suppose that a randomly selected 
subsample of 300 clusters from the whole Planck catalogue 
have been observed optically and that we know the redshift 
for each cluster in this subsample (see section ^ . The input 
model used to simulate both data sets was the ACDM model 
in table 1. That model is compatible with the mass function 
given in Bahcall & Cen (1993), the temperature function of 
Henry & Arnaud (1991), the luminosity function of Ebeling 
et al. (1997), and the fiux function of Rosati et al. (1998) 
and De Grandi et al. (1999), as it was shown in Diego et al. 
(2000). 

We have compared both simulated data sets with ~ 2 mil- 
lions diff'erent flat ACDM models where the six free param- 
eters of our model have been changed on a regular grid. The 
first three parameters correspond to the cosmological ones 
((Tg, n, and r) which control the cluster population in the 
Press-Schechter formalism. The other three parameters cor- 
respond to the T — M relation (eq. ^ whose free parameters 
will enter in the fitting procedure at the same level as the 
cosmological ones. This relation is needed to build the to- 
tal flux from the mass of the cluster (eq. ^ . By considering 
the T — Af as a free parameter relation, we will check the 
influence of the uncertainty in the scaling relation on the 
determination of the cosmological parameters. 

In fig. ^ we show the results of our fit. We have marginalized 
the probability over each one of the six free parameters. The 
probability was defined as in Diego et al. (2000) using the 
Bayesian estimator given in Lahav et al. (1999). 

^2lnPL=xi, (9) 
where, 

2 

i 

Xi is the ordinary for each data set and Ni represents the 
number of data points for the data set i. Based on a Bayesian 
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Figure 8. Marginalized probability in each parameter. The black 
dot represents the input model (ACDM model in table 1). 



approach with the choice of non-informative uniform priors 
on the log, those authors have seen that this estimator is 
appropriate for the case when different data sets are com- 
bined together, as it is our case. The factor Ni plays the role 
of a weighting factor. Larger data sets are considered more 
reliable for the parameter determination. 
As shown in that figure, the estimate of the cosmological pa- 
rameters is unbiased (compare with the input model, black 
dots). They are also very well constrained with a small de- 
generacy between the parameters. This result shows how 
with future SZE data it will be possible to discriminate 
among several scenarios of cluster formation. The situation 
is different in the parameters of the T — M relation. In this 
case we do not get any spectacular result. Only the %j) param- 
eter has been well located. There is a degeneracy between 
the amplitude To and the exponent a which will be discussed 
in the next section. 

In fig. ^ and fig. ^ we present the simulated data sets and 
the two undistinguishahle models given in table 1. From 
the first figure it is evident that both models would remain 
undistinguishahle if only that data set is used in the fit but 
in the second figure we can see that it is possible to distin- 
guish the two models at a high significant level due to their 
different evolution in redshift. 



Figure 9. dN/dS curve for the ACDM (solid) and OCDM (dot- 
ted) models in table 1. Data points represent a Monte Carlo re- 
alization of mock data for the ACDM model. This curve includes 
all the Planck detected clusters (2/3 of the sky). 
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6 DISCUSSION & CONCLUSIONS 

In previous sections we have seen that the SZE will be a 
very powerful tool to study the cluster population at dif- 
ferent redshifts. Up to now, no cluster has been detected 
above z « 2.0. Previous X-ray surveys have been limitted 
in redshift and current experiments (CHANDRA, Newton- 
XMM) are not expected to detect clusters much above that. 



Figure 10. dN{S > 30mjy,z)/dz (353 GHz) curve for the 
ACDM (solid) and OCDM (dotted) models in table 1. As in fig- 
ure ^ the data was simulated assuming a ACDM model. The 
simulation was over a solid angle dQ = 0.6675% of the sky where 
a total of Si 300 clusters were found with fluxes greater than 30 
mjy. 
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Only through the SZE we can have the possibiUty to observe 
clusters above that redshift (or maybe to conclude that no 
cluster has been formed above that redshift). These high 
redshift clusters are fundamental to understand the physics 
of cluster formation and also to establish the evolution of 
the cluster scalings such as the T — M ov — M . 
We have seen that Planck will be able to detect distant clus- 
ters which will provide very useful information about the 
cluster population and the underlying cosmology. 
However, we have seen that with only the dN{S)/dS curve, 
it will be difficult to discriminate among models which were 
previously undistinguishable. To distinguish them, we need 
redshift information. We have seen that for the whole SZE 
sky catalogue, only a relatively small number (« 300) of op- 
tically observed clusters randomly selected from the whole 
Planck catalogue is needed in order to discriminate between 
the ACDM and OCDM models just by looking at the differ- 
ent cluster population as a function of redshift. If we want 
to discriminate among the ACDM models, we have shown 
that by combining the statistically large data set dN{S)/dS 
with the cosmological sensitivity of dN{> SOmJy, z)/dz it 
is possible to reduce significantly the degeneracy in the cos- 
mological parameters as can be seen in fig. (^. 

One important conclusion is that this result is almost inde- 
pendent of the assumed T — M relation. In fact our method 
is practically insensitive to the To amplitude and a exponent 
in eq. We have marginalized the probability assuming 
different fixed values for To and a. The resulting marginal- 
ized probabilities were very similar in all the cases consid- 
ered showing the small dependence on the assumed values 
of To and a. The almost null dependence on a can be easily 
understood by looking at eq.M. In the computation of the 
temperature function (see eq.^) the derivative dM/dT was 
inversely proportional to aM°'~^ . The X-ray derived func- 
tions (like the temperature function) are sensitive to the a 
exponent through the previous derivative. On the contrary, 
the flux function, d? N / dSdz is inversely proportional to the 
derivative dS/dM oc (1 + a)M°'. Therefore, a change of 0.1 
units in a, represents a change in the dM/dT derivative of 
14% while in the flirx function the same change in a implies a 
variation of only 6% in the derivative dS/dM, both percent- 
ages assuming M — 1 x 10^^ Mq. This explains why the 
SZ flux function is less sensitive to a than the X-ray derived 
functions. The uncertainty in To is a bit more difficult to un- 
derstand. From eq. (P) the total flux is directly proportional 
to To and therefore we should expect some dependency of 
our fit on this parameter. However, if a change in To is com- 
pensated by a change in a then we would have a degeneracy 
on these two parameters (S cx roM^+"(l -1- z)''' /Daizf). In 
fact from fig. 0we can see that those models with a low value 
for To and a high value of a are slightly favored indicating 
this fact that there is some kind of compensation between 
these two parameters. 

In order to break the degeneracy in To — a we should in- 
clude in our fit information concerning the mass of the clus- 
ters just to make the fit sensitive to an independent change 
in To and/or a and not only to a change in the quantity 
ToM^+". The previous situation was considered in Diego et 
al. (2001) where we included in the fit the cluster mass func- 
tion. In that case we found in fact that there was not any 
degeneracy in those parameters. 



The third parameter of the T — M relation (tp) seems to be, 
however, very relevant for our fit. This is not surprising as 
we are using one data set which is expressed as a function 
of redshift (fig. [o|). While both. To and a can be mutually 
compensated, the effect of changing on the simulated data 
sets (figs, ^and |l^) can only be compensated with a change 
in some of the cosmological parameters (through their effect 
on the cluster population and in Da{z)) but as the allowed 
range of variation of the cosmological parameters is small 
(see fig. consequently the confidence interval for ip will 
be small as well. 

In this work, the T — M relation was considered as a free 
relation just for consistency with our previous work. When 
fitting SZ data, we have shown that the choice of one spe- 
cific value for To and a in the T — A4 relation is not quite 
relevant, although it is important to include in the fit the 
possible dependency of this relation with z. This situation is 
opposite to the one in Diego et al. (2000) where the redshift 
dependence was not relevant (since most of the data was at 
low redshift) but the choice of To and a was important to 
obtain a good fit to the X-ray and optical data considered 
in that work. The specific form of the T — M relation will 
be more important in the case of fitting future X-ray data. 
Newton-XMM will provide very relevant information, spe- 
cially at low and medium redshift, about the cluster popula- 
tion and the scaling relations T—M and L^ — T. However, we 
have seen that the expected number of detected clusters and 
the redshift coverage will be smaller for this mission com- 
pared with Planck and therefore Planck will provide several 
key informations to understand cluster formation and evolu- 
tion. For instance, as we have already seen, the information 
about the T — M relation can be complemented with studies 
of the SZE on clusters. Meanwhile To and a can be deter- 
mined through the study of low redshift X-ray data, could 
be constrained with the high redshift SZE data. The best re- 
sults will come, therefore, from the combination of data from 
X-ray and mm missions (see eg. Haiman, Mohr & Holder, 
2000). With Newton-XMM we can obtain a good sampling 
of the cluster population at low-intermediate redshift with 
their corresponding temperatures and X-ray fluxes (also de- 
tecting the low mass population) and with Planck we will 
explore the cluster population further in redshift. 
A very interesting possibility has been analyzed by Xue & 
Wu (2001). The authors suggested the use of the X-ray lumi- 
nosity function as a starting point to derive the mm (SZE) 
flux function. In the process, several assumptions about the 
Lx ^ T and T — M relations need to be done. These as- 
sumptions could be tested with future SZE data opening 
the possibility to study those relations at redshifts where no 
clusters can be observed in the X-ray band. 

Although this paper has concentrated on the possibilities of 
the future Planck SZE catalogue, proposed and undergoing 
mm experiments will measure the SZE for hundreds of clus- 
ters before Planck is launched. These experiments will open 
a new era in which many works will be done based on those 
exciting data. 
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APPENDIX A: PLANCK FLUX LIMIT 

In this appendix we want to justify that Planck will be able 
to detect those clusters with an integrated total flux above 
^ 30 mjy. 

Obviously, the number of detected clusters will depend on 
the technique used to detect them. We will focus our atten- 
tion on the Maximum Entropy Method (MEM) (Hobson et 
al. 1998, Hobson et al. 1999) where the authors have shown 
that with such method they obtain a good recovery of the 
thermal SZE. In that paper it is shown that the rms resid- 
uals per 4.5 arcmin FWHM Gaussian beam for the MEM 
reconstruction is ~ 6/iA" per pixel. 

Now we compute the flux (in mJy) corresponding to that 
rms temperature and therefore it should be considered as a 
flux per pixel. 

The flux is defined as the integral of the specific intensity 
on the solid angle 



Ssze{u) = / AIszE{iy)dQ 



(Al) 



When we compute the total fiux of the cluster the solid angle 
is that subtended by the cluster (see eq. (^) in section g). 
In eq. ([ai|), AI{iy) can be related with AT/T as AJ(i/) = 
Iog{x){AT /T) where lo is a constant given below, g{x) is 
the spectral shape factor, AT — 6.0 fiK is the temperature 
we want to transform into a fiux and T — 2.73 K, 



{hp 



= 2.7 X 10' 



mJy 

str 



(A2) 



(A3) 



And the spectral dependence g(x) is given by: 

Q{x) — -, 

> (go; „ 1)2 

where x = hv/{kBTcMB) ~ uiGH z) /56.8 is the adimen- 
sional frequency. Assembling terms in eq. (Al) we get the 
(1.5 arcmin)^ = 1.9 x 10" 

mJy 



flux per pixel (dil 
Ssze{i^) 



^ str): 



AT 

2.7 X 10^^— g{x)dn- 
"^'^y (SOOGHz) 



str 



(A4) 



0.45- 



pix 



This number has been calculated for the frequency v — 300 
GHz. In the paper we presented our calculations for 353 
GHz. The flux at this frequency is a factor « 1.3 times 
higher than at 300 GHz. Therefore 

mJy 



SszE = 0.58- 



pix 



{353GHz 



(A5) 



This fiux should be considered as the rms in the residual map 
when the SZE is recovered by MEM, that is the noise per 
pixel. We will consider a conservative limit for the detection 
of a cluster of signal to noise ratio > 7.5a on the FWHM. If 
we consider our antenna as a Gaussian beam and the cluster 
is not resolved (this will happen in most of the cases in 
Planck) then we can consider that the cluster profile follows 
a Gaussian pattern: 



Signal = S = Ac 



V2a 



(A6) 



Inside the FWHM the signal is just the integral of the pre- 
vious equation from r = to r = FWHM/2 which can be 
solved easily (by doing the variable change x — r^ /2a'^): 
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(•FWHM/2 

S{FWHM) = A I = 6.3 x A (A7) 







Inside the FWHM the noise associated to MEM goes hke: 



N{FWHM) = Nr^s \J Nl'^„^,j ( A8) 

where Nrms is the noise per pixel found previously and 
^FWHM ^® number of pixels corresponding to the area 
enclosed by the FWHM, iVl^'^^jv^ = FWHM^ = 11.1 
(FWHM = 5 arcmin for the channel at 353 GHz). 
Now if we require S(FW HM) /N(FW HM) = 7.5 then 



by dividing eqs. (A7) by (A8) we get that A must be 



A ~ 4 X Nrm s, th at is, an unresolved cluster with a signal 



following eq. (A6) should have an amplitude A 4 x Nrms 
in order to have S/N = 7.5 inside the FWHM. 
The total flux inside the full antenna beam of such a sig- 
nal is Srotai ~ 12.6 X j4 « 50 X Nrms- If now we substi- 
tute Nrms ~ 0.58mJy per pixel, then we finally obtain, 
Srotai ~ 29mJy which approximates the value of 30 mjy 
used in the paper. 

This paper has been produced using the Royal Astronomical 
Society/Blackwell Science ET]eX style file. 
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